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ABSTRACT 

A simple analytical model is used to derive the main properties of super- 
novae type la (SNe la), which are produced by the thermonuclear explosion 
of accreting C-0 white dwarfs that cross the Chandrasekhar mass limit. The 
few underlying physical assumptions of the model yield analytical expressions 
that reproduce quite well the observed bolometric light-curves of SNe la and the 
empirical brighter-slower and brighter-bluer relationships that were used to stan- 
dardize SNe la and use them as distance indicators which led to the discovery of 
the accelerating expansion of the universe. 

Subject headings: supernovae: general 



1. Introduction 

Despite large observational, theoretical and numerical efforts over decades, supernovae 
explosions are not fully understood. Standard stellar evolution theory predicts that stars 
with initial mass M^8 Mq end their short life in core-collapse supernovae explosions, while 
long-lived stars of mass less than ~ 8 Mq eject most of their outer layers, leaving a ~ 1 Mq 
carbon oxygen (C-0) white dwarf (WD) that cools slowly by radiation, possibly for billions 
of years. If such a C-0 white dwarf accretes mass from a companion star and crosses the 
Chandrasekhar mass limit, Mch ~ 1.38Mq (Chandrasekhar 1931), at which the pressure of 
its degenerate electron gas can no longer balance the gravitational pressure, it undergoes a 
gravitational collapse. Its temperature rises sharply and triggers a runaway thermonuclear 
explosion (Hoyle and Fowler 1960) where no central object is left over, most of the nuclear 
binding energy release is converted to kinetic energy of the debris of a supernova (SN) type 
la, while its observed bolometric light-curve is powered mainly by the decay of the long-lived 
end-product radioactive elements that were synthesized in the explosion. 

As the progenitors of core-collapse SNe are much more luminous than WDs, they have 
been identified from archival images on multiple occasions Smartt 2009, while there have 



^Physics Department, Technion, Haifa 32000, Israel 



- 2 - 



not yet been any direct observations of the WD progenitor of an SN la. However, the above 
paradigm of the origin of SNe la is strongly supported by several observational facts (see, 
e.g., Astier 2012 and references therein) such as: 

• SNe la are the only type of SNe observed in old stellar environments such as eUiptical 
galaxies. WD are produced by long-lived stars of less than ~ 8 Mq that eject most 
of their outer layers, leaving a ~ 1 Mq C-0 WD that cools slowly by radiation, for 
billions of years. 

• The fast exothermic nuclear fusion reactions, starting with helium capture by ^^C and 
^^O, produce the intermediate mass elements Si, S and Ca observed in the spectrum 
of SNe la but not H and He, which are lacking in the initial state and in the spectra 
of SNe la. 

• The fast exothermic nuclear fusion reactions end with the near center production of ^^Ni 
whose radioactive decay chain ^^Ni(r = 8.76 d) Co(t = 111.27 d) -^^^ Fe seems to 
power the light-curves and explain the late appearance of iron group elements' spectral 
lines. 

• The kinetic energy of the debris, the bolometric light-curve, the spectrum and the 
spectral evolution of SNe la are roughly those expected. 

Many impressive hydrodynamical codes with radiation transport that calculate the light- 
curves and spectra of SNe have been developed in the past two decades. Such codes require as 
an input the initial distribution of density, chemical composition, velocity, temperature and 
angular momentum in the SN ejecta. These, in principle should have been provided by the 
numerical simulation of the initial explosion. However, so far such explosion codes starting 
with configurations obtained from standard stellar evolution calculations have not been able 
to produce consistently supcrnovae explosions of any type. Nevertheless, it is widely believed 
that such elaborate codes are necessary in order to reproduce the light-curves and spectra of 
SNe of all types during their photospheric and early nebular phases (see, e.g., Frey et al. 2012 
and references therein). This is because of the enormous complexity of SN explosions and 
the large diversity of their light-curves, complex spectra and spectral evolution. 

Despite the complexity of SN explosions, SNe la light-curves and spectra display an 
approximate 'standard candle' behaviour (see, e.g.. Branch & Tammannl992) with simple 
correlations between various photometric properties such as peak luminosity, color and de- 
chne rate following peak luminosity (e.g., PhiUips 1993; Riess et al. 1996; Tripp 1998). These 
correlations were used to standardize SNe la and have allowed to improve the precision of 
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cosmic distances estimated from SNe la observations that led to the discovery of the accel- 
erating expansion of the universe (Riess et al. 1998; Perlmutter et al. 1999). The nearly 
standard properties of SNe la, however, suggest that perhaps the bolometric hght-curves 
and other general properties of SNe la can be derived directly from a simple model, despite 
the complexity of SNe la explosions and the details of the complex radiation transport in 
the ejecta. Semi-analytic approaches have been pioneered by Colgate and McKee (1969), 
Arnett (1979), and Colgate, Petschek and Kricsc (1980). Improved analytic solutions for the 
conversion of radioactive decay energy into the light-curves of SNe la have been proposed 
more recently, e.g., by Pinto and Eastman (2000). Here, using an even simpler analytical 
model, we derive the main standard properties of SNe la during their photospheric phase. In 
particular, we derive simple closed form expressions that well describe their bolometric and 
photometric light-curves during their photospheric and early nebular phases. In a following 
paper we extend that to the main types of core-collapse SNe, SNe IIP, SNe IIL, SNe Ib/c, 
Super Luminous SNe of types 1 and 11 (see, e.g.. Gal- Yam 2012 and references therein for 
the current classification of SNe), which are powered mainly by the radioactive decay of ^^Ni 
and/or by the interaction of their ejecta with their environments. 



2. Derivation of SNe la General Properties 

We adopt the current paradigm of SNe la explosions, where a carbon-oxygen (C-0) white 
dwarf accretes mass from a companion until it reaches a critical mass Mf. {M^ ~ 1.38 Mq, 
the Chandrasekhar mass limit) and undergoes a gravitational collapse that rapidly raises its 
temperature and triggers thermonuclear reactions that produce intermediate mass elements 
(IME=Si,Ca,S,...) via deflagration (Nomoto et al. 1976, 1984) followed by a transition 
to detonation that converts a fraction rj of them to ^^Ni whose decay chain powers the 
bolometric light-curves of SNe la (Colgate & McKeel969; Colgate et al. 1980; Kuchner et 
al. 1994). This fraction 77 probably depends on the unknown initial conditions (such as 
distributions of mass, composition and angular momentum in the WD progenitor when its 
mass crosses M^), and will be treated as a free adjustable parameter. We assume spherical 
symmetry and ignore environmental effects produced, e.g., by a companion star or a common 
envelope. 

(a) The kinetic energy of the explosion: The nuclear binding energy release per nucleon 
in the explosion when a C/0 nucleus is synthesized into IME is f« 8.551 — 7.828 = 0.723 MeV 
or ?a 8.642 - 7.828 = 0.814 MeV if the C-0 is synthesized into ^"^Ni. A negligible fraction 
of these energies escapes through photons and neutrinos and most of it is converted to the 
kinetic energy of the explosion. Morover, when the mass of an accreting WD crosses Mc, 
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the sum of the gravitational binding energy and free energy of the degenerate electron gas 
is 0, independent of radius. Consequently, the kinetic energy Ej^ of the explosion at oo is 
approximately the nuclear binding energy released in the synthesis of IME and ^^Ni, 

Ek « Nb [0.723 + f] 0.091] MeV (1) 

where is the baryon number of the C-0 WD progenitor. For Mc — Men — 1.38 M©, 
Nb ^ 1.67 X 10^^ and E^ ^ (2.04 ± 0.12) x 10^^ ergs. 

(b) The expansion velocity: We neglect angular momentum and enviromental effects 
and assume an homologous expansion of the SN fireball, i.e., a uniform spatial density 
throughout the fireball at any moment. This implies that the expansion velocity v{r) at 
distance r from the center satisfies v{r) — (r/R) Vexp where R is the radius of the fireball 
and Vexp is its expansion rate at R. Consequently, the total kinetic energy of the explosion 
is Ek = (3/10) NbUibV^^p. Substituting that into Eq. (1) yields 

Veccp = VlO (0.723 + ri 0.091) MeV/3 nib (2) 

where the effective mass of a baryon in ^°Ca and ^^Ni is ^ 930.48 MeV/c^. Hence, the 
expected initial expansion velocity is roughly Vg^p ~ 15, 740 ± 460 km/s independent of the 
value of Mc- Because of the typical low density of the surrounding medium where SNc la 
explosions take place, the deceleration of the expansion due to swept-in mass is negligible 
during the photospheric and early nebular phases, unless the explosion takes place inside a 
common envelope of a close binary. Note that as long as the SN fireball is opaque to its black 
body radiation, the observed photospheric velocity Vph is roughly the expansion velocity at 
a depth where the optical opacity satisfies T{Rph) ~ 2/3, i.e., Vph ^ {1 — 2/3t)T4xp, where 
T is the optical depth to the center of the fireball. 

(c) The bolometric luminosity: Let t — he the time of shock break-out. The typical 
environment of SNe la has too low a density to slow down the expansion of the SN fireball. 
Consequently, for t > the radius R of the SN fireball expands homologously with practically 
a constant R — V^^p- The bolometric light-curves of SNe la are powered by trapping the 
energy of the gamma rays and positrons from the radioactive decay chain ^^Ni -^"'''^ Co — ^►^^ Fe 
near the center. Throughout the photospheric phase, the SN fireball is highly opaque to 
radiation and the energy input is mainly through Compton scattering of the gamma-rays. 
The thermal energy density u{T) is dominated by black body radiation, i.e., u{T) aT^, 
where a — 7.56 x 10~^^ergcm~^K~^, and the total thermal energy of the SN fireball is 
U ^ {An R^/3) aT^. During the photospheric phase, the SN fireball loses energy mainly by 
expansion, at a rate ~ PdV/dt U/t (for a constant K^p); a-nd by emission of photons, 
which are transported to the surface by a random walk. The photon emission yields a 
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bolometric luminosity L ^ U/tdif, where tdif ~ R^/Xc — Rt/c is the mean diffusion time 
of photons to the surface by a random walk, A is their mean free path and t — RuiUi is the 
fireball opacity where the summation extends over all particles in the fireball (ions, neutral 
atoms and free electrons) with density rii and effective cross section (Tj. Roughly ~ R~^ 
and R = Vexpt. Hence, the mean diffusion time decreases with time roughly like t^^ and 
can be written as tdij ~ t^/t where tr oc mI^"^ . During the photospheric phase, energy 
conservation can be approximated by 

U + U[\ + ^]^Era, (3) 

where E^a is the energy deposition rate in the fireball after the thermonuclear explosion by 
the decay of radioactive elements synthesized in the explosion. The solution of Eq. (3) is 

U^ - / te'"'^'"Eradt. (4) 

t Jo 

Consequently, the bolometric luminosity that satisfies Lb — tU/t'^ is given by the simple 
analytic expression 



f te'^'^'^Eradt. (5) 



For a luminosity that is powered by the radioactive decay chain ^^Ni — >-^^ Co — >-^^ Fe, 
Era — E^ + Ee+ where 

E^ = M!!^ [A^^Ni) 7.85x10^^ e-'/^-'^'+A^iCo) 1.43x10^^ [g-t/iii.27d_g-t/8.76d] ^^g^-i _ 

Mq 

(6) 

Aj{Ni) and Aj{Co) are the attenuation of the 7-rays from the decay of ^^Ni and ^^Co, 
respectively, in the SN fireball. A^^ ^ 1 — e~'^'' , where r^(t) ~ i? < njO",; > is the fireball 
opacity and the summation extends over all its particles. For a composition dominated by 
IME and iron group elements, pi Aj ~ 12gcm^^ for the 7 rays with < E'^ 0.54 MeV from 
the decay of ^^Co and larger by a factor 5/3 for 7-rays with < E^ >^ 1.32 MeV from the 
decay of ^^Ni (see Hubbell 1982). Hence, = Rp/pX « tyt^ where « 3 MJ4: tt V^^ pi A^, 
and t^{Ni) ?a 55 d while t^(Co) ^ 42 d, for ^ Men, 

The positrons from the /3+-decay of ^^Co (and the e"^ from the decay of other relatively 
long lived radioactive isotopes that were synthesized in the thermonuclear explosion) are 
presumably trapped by the turbulent magnetic field of the SN fireball. The power supply 
by positrons and e'^e~ annihilation into two 0.511 MeV 7-rays is approximated by 

E,i = (1 + 6.22) A^) 1.43 x 10^^ [e-^n^^-^^^ - e-*/8.76d] ^^g^-i _ 

Mq 
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where Ag is the ratio between the energy released in the radioactive decay of ^^Co as kinetic 
energy of e'^ and that released directly in 7-rays. The kinetic energy deposition dominates the 
power supply when the fireball becomes highly transparent to 7-rays, i.e., when t ^ t-y{Co). 
During this nebular phase the ionization and excitations by 7-rays and lead to scintillation 
and bremsstrahlung emission that dominate the SN spectrum. 

(d) Ecirly- and late-time behaviours of Lb{t): Era changes rather slowly with t relative 
to the fast rise of te* Z^*'- with increasing t, and can be factored out of the integration in 
Eq. (5) yielding ^ [1 — e~*^/^*^] ij^-a. Hence, the bolometric luminosity rises initially 
like Lij ^ (t^ /2t1) Era and has the late-time asymptotic behaviour Li, — > Era- Note that 
the derivation of Eq. (5) is valid only for the photospheric phase. However its late-time 
behaviour Li-,(t) ^ Era(t) is also the correct behaviour of Lf, during the nebular phase. Thus 
Eq (5) is valid for both the photospheric and nebular phases. 

(e) The colour temperature during the photospheric phase: As long as the fireball is 
optically thick (r ^ 1), its continuum spectrum is approximately that of a black body, i.e., 
I,y (X j[e^^/^'^ — 1), where the effective temperature T satisfies the Stefan-Boltzmann law. 
For homologous expansion, the photospheric velocity decreases like Vph ~ ~ 2/3r)), 
and the Stefan-Boltzmann law yields an effective temperature 



T 



1/4 

(8) 



where a = 5.67 x IQ-^ ergs"^ cm'^ K^^ Thus, T ^ [^^„(t)/4 tt al/y t^^/^ fo^ ^ ^2 
while T [Era{t) /'i'n' aVp^Y^"^ t~^/'^ for » (as long as r » 1). During the transition 
from the photospheric phase to the nebular phase, when free-free emission and scintillations 
take over, the temperature remains roughly constant. 

(f) Luminosity peak-time: A very rough but very simple estimate of tp follows from 
approximating the early-time bolometric luminosity by Lt,{t) oc Eraif) where E^ait) is 
dominated by the decay of ^^Ni, i.e., Era{t) oc e"*/^-'''^''. This approximate form of the early- 
time Li,{t) has a maximum at i = 2 x 8.76 = 17.52 d, which depends neither on Mc nor on 
the synthesized mass of ^^Ni. 

The approximate expression ~ [1 — e"*^/^*^] Era peaks at t = f« 17 ± 1.5 d for 
ij. 11 ± 1 d. This peak-time does not deviate significantly from that of Li,{t) as given by 
Eq. (5). It depends on i^, i-e., on since tr oc M^^^, but not directly on the synthesized 
mass of ^^Ni. 

(g) VeaAa luminosity - nickel mass relation: At its peak-time, t — tp, the bolometric 
luminosity Li,{t) = t U/tl satisfies = U/tl + ttf /tr ^ 0. It then follows from Eq. (3) that 
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Lb{tp) = Era{tp). At t — tp the SN is still highly opaque to the 7-rays and hence 

L,{tp) ^ ^^2^'^ 6.35 X 10^^ [e-*-/«-^6'^ + 0.236 e-*-/iii-2^'^] ergs'S (9) 

which yields Lb{tp) ^ (2.2 ± 0.010) x 10^^ (^M{^^Ni) j Mq) ergs"^ for = 17 ± 1 d. 

(h) Optical and UV light-curves and peak times: As long as the SN fireball is optically 
thick, it radiates like Planck's black body and the light-curves at a frequency v satisfy 

where h is the Planck constant and k is the Boltzmann constant. Near the peak-time of the 

bolomctric luminosity, the temperature has the approximate behaviour T(t) ~ T{tp) [tpjtY^'^ ■ 
Morover, e^^l^^^ 1 in the VBU bands and since Rph act ior t <^ t^, the maximum of L^{t) 
is reached there aX, t — tp^i, where 

tp,^ ^ WikTp/huftp oc v^MpTVi) . (11) 

(i) The decline rate measure Ami5(A): Although the peak intrinsic luminosity of SNe 
la is not a standard candle, it appears to be correlated to the shape of their light-curves 
(Phillips 1993). Since 1993, various empirical correlations between the peak absolute magni- 
tude of SNe la and the measured shapes of their intrinsic light-curves have been adopted for 
the use of SNe la as standard candles for distance measurements (e.g., Hamuy et al. 1996; 
Riess et al. 1996, 1998; Perlmutter et al. 1999; PhiUips et al. 1999; Goldhaber et al. 2001; 
Prieto et al. 2006). Most methods have used Ami5(A), the magnitude difference between 
maximum brightness and 15 days past it, as a measure of the decline rate. For e^^^^'^'p 3> 1 
and ~ 11 d (ip « 18 d), it follows from Eq. (10) that 

he 

Ami5(A) -1.28 + 0.52— —. (12) 

\ K Ip 

Because Tp oc [M(^^iVi)]^/'^, Ami5(A) decreases with increasing Mi^^Ni) and produces the 
'brighter-slower' behaviour of colour hght-curves. Note also that the SN fireball during 
the photospheric phase is essentially Planck's black body whose radiation satisfies Planck's 
radiation law (Eq. 10), it has the usual 'brighter-bluer' behaviour of a black body radiation. 



3. Compcirison with observations and conclusions 



The properties (a)-(i) of standard SNe la that were derived here from a simple model are 
in very good agreement with those observed. This is demonstrated in particular in figures 1-3, 
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where we compare our analytic expression Eq. (5) and the observed rest-frame bolomwtric 
hght-curves of three nearby SNe la with very early detection and continuous follow up: 
SN2011fe in MlOl (Nugent et al. 2011; Munari et al. 2012), SN1994ae (Contardo et al. 2000) 
and SN1992bc (Contardo et al. 2000). The spread in arrival times At < R/c ^ Vpht/c < t/20 
of simultaneously emitted photospheric photons was neglected. The best fit values of the 
4 adjustable parameters are listed in Table I. Their values are consistent with their simple 
estimates in the text. The empirical brighter-slower and brighter-bluer behaviours of SNe la 
light-curves that allowed their use as standard candles for distance measurements and led 
to the discovery of the accelerating expansion of the universe, are satisfied by the predicted 
light-curves. 

SNe la, like all other types of SNe, are very complex phenomena. They depend on 
the mass of their WD progenitor, on the detailed 3D distributions of its density, chemical 
composition and angular momentum and on its environment, which are largely unknown. 
These unknowns determine in a complex way the hydrodynamics and radiation transport 
in SN explosions. Nevertheless, here we have demonstrated that the main properties of 
SNe la can be well predicted from simple considerations and their lightcurves during their 
photospheric phase can be well described by simple analytic expressions which involve only 4 
adjustable parameters {M{^^Ni), tr, t^, Af,), whose values do not deviate significantly from 
those expected on general grounds. 
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Table 1. Best fit parameters of tlie analytic description, Eq. (5), of the bolometric 
light-curves of the three SNe la shown in Figs. 1-3. to is the explosion time relative to 

maximum light 



SN to [d] U [d] t^{m) [d] M(56Ni) 



SN2011fe -17.07 13.16 24.05 0.732 Mq 0.033 

SN1994ae -16.03 9.62 37.64 0.454 0.042 

SN1992bc -18.66 11.29 35.31 0.782 Mr,, 0.040 
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Fig. 1. — Comparison between the bolometric light-curve of SN2011fe (Munari et al. 2012) 
and that predicted by the analytic model and summarized in Eqs. (5)- (7), assuming it was 
powered by the decay of ^^Ni, synthesized in the thermonuclear explosion of a C-0 white 
dwarf of a critical mass. 
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Fig. 2. — Comparison between the measured bolometric light-curve of SN1994ae (Contardo 
et al. 2000) and that predicted by the analytic model and summarized in Eqs. (5)- (7), 
assuming it was powered by the decay of ^^Ni, synthesized in the thermonuclear explosion 
of a C-0 white dwarf of a critical mass. 
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Fig. 3. — Comparison between the bolometric light-curve of SN1992bc (Contardo et al. 2000) 
and that predicted by the analytic model and summarized in Eqs. (5)- (7), assuming it was 
powered by the decay of ^^Ni, synthesized in the thermonuclear explosion of a C-0 white 
dwarf of a critical mass. 



